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Abstract. We present new observational constraints on inliomogeneous models based on ob- 
servables independent of the CMB and large-scale structure. Using Bayesian evidence we 
find very strong evidence for the homogeneous LCDM model, thus disfavouring inhomoge- 
neous models. Our new constraints are based on quantities independent of the growth of 
perturbations and rely on cosmic clocks based on atomic physics and on the local density of 
matter. 
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1 Introduction 

The discovery of accelerated expansion about fourteen years ago, based on the dimming of 
distant supernovae (SN) ([1, 2]), has led to a standard model of cosmology (ACDM) in which 
about 72% of the energy of the universe is in the form of a cosmological constant (CC) (or 
dynamical dark energy). An important assumption in this model is that our universe is 
homogeneous. 

As an alternative to a cosmological constant, which has significant theoretical problems 
due to the required fine-tuning, the supernova data could also be explained without a CC 
if we live close to the center of a large (> 1 Gpc), spherically symmetric void (see [3-5] for 
some of the first such proposals, while [6, 7] discussed this idea even before the measurement 
of accelerated expansion). The reason is that we cannot distinguish the effect of spatial 
variations in the geometry from temporal ones in the radial direction from measurements 
along the past light-cone alone. Specifically, in a large void, the local expansion rate gets 
larger closer to the center of the underdensity, thus causing the same additional dimming at 
large distances that would be caused by accelerated expansion in a homogeneous cosmology. 

Such void models arguably have larger philosophical and aesthetic problems than the 
CC. First, the existence of a void of the required size is extremely unlikely in the standard 
inflationary scenario and its generation would probably involve rather unusual early universe 
physics. Moreover, there is a large amount of fine tuning associated with the requirement 
from observation that we find ourselves within ~ 1% of the void's center ([8-13]), a grave 
violation of the Copernican principle. Additional fine-tuning comes from the fact that the 
void has to be close to spherical to be consistent with the isotropy that we observe. Finally, 
the void as an explanation for cosmic acceleration does not actually solve the cosmological 
constant problem. It merely assumes that the CC is equal to zero without providing any 
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explanation (this criticism of course also applies to most if not all models of dynamical dark 
energy) . 

Still, the final say can only come from observation so it is worth comparing predictions 
of the void model to recent cosmological data. This is not just interesting for the sake of 
understanding the viability of these specific models, but also serves to test a not often directly 
tested pillar of our cosmological framework, the assumption of homogeneity (see e.g. [14]). 

The void is commonly modeled as a Lemaitre-Tolman-Bondi universe (LTB, [15-17]), 
which describes the most general spherically symmetric universe filled with pressureless mat- 
ter (dust). Such models have been extensively used to reproduce the observed supernova 
magnitudes and can in fact explain any Hubble diagram perfectly, provided enough freedom 
in the void profile is allowed (see e.g. [7, 18]). A more trying test of the LTB cosmology arises 
when additional data sets are added. Many data sets have by now been considered, including, 
but not limited to, the primary cosmic microwave background temperature power spectrum 
(e.g. [19-22]), spectral distortions of the cosmic microwave background (CMB) ([23]), the 
kinetic Sunyaev-Zel'dovich (kSZ) effect ([24-28]), the primordial Lithium abundance [29], 
the Baryon Acoustic Oscillation (BAO) scale (e.g. [20, 30]), and combinations of the above 
([10, 20, 21, 31-33]). In particular, both the combination of the CMB with a low redshift 
Hubble parameter measurement (and with other data sets) ([31]), and the observational up- 
per limit on the strength of the kinetic SZ effect in the CMB ([24-28]), appear to rule out at 
least the simplest versions of the void cosmology. 

However, many of the previous studies rely either on the CMB or, through the BAO, 
on large scale structure (LSS), while it is not universally accepted that these observables are 
understood well enough in an inhomogeneous universe. Particularly, solving the evolution 
of density perturbations in LTB is notoriously difficult, although significant progress has by 
now been made ([34-39]). Moreover, constraints derived from the CMB power spectrum may 
depend on implicit assumptions regarding the distribution of radiation ([29, 40]). 

It is thus worthwhile considering constraints that do not depend on either CMB or LSS, 
nor on the details of perturbation evolution in general. For this reason, we will constrain 
the geometry /expansion of the universe with data complementary to supernovae, by using 
"red envelope" galaxy ages ([41, 42]) as lower bounds on the age of the universe at various 
redshifts in the range z = — 1.9. These estimates of oldest galaxy ages in samples of passively 
evolving, massive, red galaxies have been used previously to measure the Hubble parameter 
H{z) as a function of redshift (see [41-44]). In fact, this measurement of H(z) has even 
been used to constrain LTB models in [45, 46]. However, the Hubble determination crucially 
depends on the assumption that the average formation time of the oldest galaxies in each 
sample is homogeneous. Since this is not guaranteed in a void model (see Appendix A), an 
issue which was at least discussed in [46], we cannot safely estimate H(z) from differential 
age measurements, and use the ages themselves instead. As it turns out, the supernova 
magnitudes, combined with a local determination of the Hubble rate, prefer a void in which 
the age of the universe as a function of redshift is low compared to these galaxy ages (see also, 
e.g., [33]) so that the addition of the age estimates places strong constraints on the model. 

We thus carry out a joint analysis of these three data sets, considering a void that is 
asymptotically flat and with homogeneous big bang time. We derive parameter constraints 
and, performing Bayesian model selection, find that ACDM is significantly favored over LTB 
models. The data prefer a very low value of the relative matter density in the center of the 
void (r^in < 0.1), and we show that adding a direct measurement of from, e.g., clusters 
causes additional tension between the data sets, helping to further disfavor the LTB model. 
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We thus find results consistent with the studies based on CMB and/or LSS, but using an 
independent method. 

We will briefly review the LTB model and introduce our parametrization in section 2, and 
describe the four types of data used in section 3. We then derive parameter constraints, and 
compare the model to ACDM using both the statistic and Bayesian evidence comparison, 
in section 4. Finally, we summarize and discuss our results in section 5. The dangers of using 
red envelope ages to estimate the Hubble parameter in void models are discussed in Appendix 
A. Throughout this article, we will use notation similar to that in [33] (and previous papers 
by these authors). 



2 The Model 



The most general spherically symmetric space-time for a universe filled with pressureless 
matter {p = 0) is given by the LTB metric, 

ds^ = dt" - ^^^^ dr^ - A\r, t) dn\ (2.1) 
1 — k[r) 

where (Kl^ = dO'^ + sin^ 6d(jP' is the squared line element in the transverse direction. Here, 
k{r) is the (radial) position dependent curvature function and A{r, t)/r a generalized, position 
dependent scale factor. Throughout this work, a ' denotes a partial derivative with respect 
to radial coordinate r and a ' will denote a partial time derivative. The LTB metric allows 
for different rates of expansion in the longitudinal and transverse directions, quantified by 
Hl = A! I A! (the longitudinal Hubble rate) and Kq^ = Aj A (the transverse Hubble rate), 
which are not necessarily equal. 

The Einstein Equations imply a generalized Friedman equation (see e.g. [30, 47, 48]), 

Hi . HSir) (..„M (^)' + ,1 - (^)^) , (2.2) 

where Hq^v) = Hx^r, to) and AQ{r) = A(r,tQ), and to is the cosmic time at present. Once the 
gauge freedom in A{r, t) related to rescalings of r is fixed, the model is completely specified 
by the two free functions Ho^r) and i^mi^) which are related to the physical matter density 
Pm{i^) and the curvature by 

J_ JdAA'pr^jr) 

^^^H,Um{r)- , (2.3) 

i.e. the average matter density within r, and by 

k{r)=Hi{r){n^{r)-l)Al{r). (2.4) 

For a given value of r, one recognizes Eq. (2.2) as the Friedman equation for a curved, 
dust-filled universe with Hubble parameter Ho(r) and matter density relative to critical den- 
sity ^Irni''')- For il.mi''') < 1, the solution is 

-^(-'-) = 2(l'^7llV)) (-M^)-^)-^o(r) 
t{r,u) - tBB{r) = H,\r) !^^^J^^^s/2 (^^^^(n) - u) , (2.5) 
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where tBB('') is in general a free function describing the location dependent Big Bang time. 
However, we will require the big bang to be homogeneous, setting tBB('^) to a constant, which 
without loss of generality we choose to be zero, 

iBB(r) = 0. (2.6) 

This condition ensures that the void can be treated as a perturbation of a homogeneous 
universe, in the sense that the void vanishes at early times ([35, 49]). 

The requirement of a homogeneous big bang restricts the freedom in model space from 
two free functions HQ{r) and 0,ra{r) to one free function and one parameter, which we choose to 
be Qmi^) and Hqq = HQ{r = 0). Once Om(?") and Hqq have been specified, to can be calculated 
at r = from Eq. (2.5) by first finding Uo{r = 0) such that A(0, Uo{r = 0)) = Aq{0) and then 
inserting uo{r = 0) into the expression for t{r,u). At any other r, requiring t{r,uo{r)) = to 
then determines Ho^r). 

In order to calculate observables in the following sections, we will calculate quantities 
along the past lightcone. The relation between t and r along a radial null geodesic is obtained 
by setting the line element (2.1) to zero, 

dt ^ ^ A'ir,t) 
dr ,Jl-k{r)' 

The photon gravitational redshift in an LTB universe is given by ([17, 47]) 

^^^ = -ff.(M). (2^8) 
We will calculate r, t and any observables of interest as a function of redshift by integrating 

-HE\r,t) 



dln(l + z) 



- ± V^lHM, (2.9) 



d\n{l + z) A'{r,t) 

and using Eq. (2.5) to find A(r,t) and its derivatives. 

Fixing Ao{r) = r, we parametrize the void profile ^Imir) by an asymptotically flat, 
constrained (in the sense that the void is not compensated by an overdensity at large radius) 
GBH profile ([30]), 

n„(,-) = l-(l-»,„)(^-''°'"l''',-f/^^^IV (2.10) 

^ ^ ^ ^ ' 1 + tanh [i2/2Ai?] J ^ ' 



The void is then completely described by the parameters 17;^ (the relative matter density in 
the center of the void), R (the void radius), AR (the width of the transition from void to 
background) and Hqq (the expansion rate at the void center). 

We will assume our galaxy is at the center of this void. If it were not, the universe 
would be anisotropic relative to our position. In particular, the observed cosmic microwave 
background dipole moment restricts us to be within about 1% of the void radius from the 
center ([8-13], but see [29, 40] for a potential caveat related to the non-negligible effects of 
radiation). 
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We have ignored the contribution from radiation in the above. Since we only consider 
low redshift observables, at times where the radiation component would be negligible, we 
expect this to be a good approximation. 

We will compare the void model described above to a homogeneous, spatially flat ACDM 
universe, filled with pressureless matter and a cosmological constant. Thus, the number of 
ACDM parameters relevant for the observables considered in this work is two and Hq), 
while in the void model, we vary four parameters. 

3 Data 

It is well known that void models can be fine-tuned to provide a good fit (in comparison with 
ACDM) to supernova luminosities. In addition, there is enough remaining freedom to provide 
a good fit to a low redshift determination of the Hubble parameter. We will therefore use 
the Union 2.1 compilation of supernova data [50] and the measurement of the local expansion 
rate by [51] as our base set of data to compare the void model against. Tension between the 
void model and the data arises when more data sets are added. The main new ingredient 
of this work is the realization that lower bounds on the age of the universe at different 
redshifts, as obtained from red envelope galaxies (see [41-44] and references therein), have 
strong constraining power. In addition, we will consider the effect of a measurement of the 
local matter density. We specifically avoid using measurements of quantities that depend on 
the details of cosmic perturbations. We describe the four data sets/compilations and their 
likelihood functions in more detail below. In the next section, we will compare these data 
against the void model and ACDM predictions, using the independence of the different data 
sets to construct the total likelihood as the product of the individual likelihoods. 

3.1 Type I A supernovae 

We use the recent Union 2.1 compilation of Type lA supernovae ([50]), which contains 557 
supernova magnitudes in the redshift range z = 0.015 — 1.4. The expectation value of the 
observed (stretch-calibrated) magnitude at the peak of the lightcurve can be written (we use 
hats to identify observables/estimators) 

{rrii) = M-51ogio {c~^H{0,to)[Mpc-^]) + 25 + blog^o [c~^ H{Q,tQ)dL{zi)) 

= M + fiiie), (3.1) 

where we have defined the distance modulus 

^ii{e) = 5 logio {c-'H{0, to)dL{z,)) , (3.2) 

which depends on cosmological parameters (and model) 9 (we have suppressed the parameter 
dependence in some places to avoid clutter). Above, M is the absolute supernova magnitude, 
which, after calibration, is assumed to be universal, 

dL{z^) = {l + zfA{r{z),t{z)) (3.3) 

is the luminosity distance to the redshift Zi of the i-th supernova and H{0, to) is the Hubble 
rate at z = and r = 0, and is equal to Hq {Hqq) in the homogeneous (inhomogeneous) case. 

Assuming the observed magnitudes rhi are drawn from a Gaussian distribution with 
means given by Eq. (3.1) and with covariance matrix Cij (we use the matrix with systematic 
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errors provided in the Union 2.1 compilation), marginalization over the nuisance parameter 
Ai (see for instance Appendix F of [52]) leads to 



(3.4) 



The normalization is such that (x^) = Nst^. Because of overlap between the Union sample 
with the data set used for the Hubble parameter described in the next subsection, we omit all 
supernovae at z < 0.1 from the Union sample to avoid double-counting of information. After 
this cut, Nsti = 392 supernovae remain. 

3.2 The local expansion rate 

We use the measurement of the Hubble parameter Hq by [51] who used Hubble Space Tele- 
scope (HST) data. This work makes use of Cepheids to construct a distance ladder to eight 
nearby Type lA supernovae, allowing a determination of the absolute supernova magnitude. 
Combining this with information from a Hubble diagram of supernovae at z < 0.1, Hq can 
be calculated, leading to Hq = 73.8 km Mpc~^ with an error = 2.4 km Mpc~^. 

This determination of the expansion rate assumes a cosmology to calculate the distance 
moduli in the supernova Hubble diagram. In ACDM, the potential error in Hq due to assuming 
an incorrect cosmology is small because at low redshift the distance moduli do not have 
a strong cosmology dependence. However, since we study cosmology beyond ACDM, we 
consider it more prudent to treat the HST measurement as a measurement of the (inverse) 
luminosity distance to the effective redshift of the supernova sample used for the Hubble 
diagram (see also [53]). The measurement of Hq given above corresponds to a measurement 



of dl^izes) = 5.97 • 10-3 Mpc"^ with a^-i = 1.9 • 10"^ Mpc'S where Zes = 0.04. The is 



Age measurements of massive, red galaxies can be used to estimate the upper edge of the 
age distribution at each redshift, the so-called red envelope ages. These measurements of the 
oldest galaxy ages vs. redshift can be used as a redshift-dependent lower bound on the age of 
the universe. In total we use 32 such age estimates in the redshift range z = 0.1 — 1.85, with 
independent error bars at the 10% level (see Fig. 3). We refer to [41, 42, 44] and references 
therein for details on the data sets used and on the age estimation from galaxy spectra. 

The likelihood of a set of (envelope) age measurements (with uncertainty o"j), given a 
cosmological model 9, is given by 



where we have not written the normalization factor, as it is independent of the cosmology 
and therefore will not be relevant for parameter constraints and model comparison. The 
label i runs over all 32 redshift in our data set, ti{6) is the cosmic time at redshift Zi in the 
cosmology 9, and j is the formation time of the red envelope galaxy at each position. In 




(3.5) 



3.3 Galaxy ages 



p{{ai}\9, {t/,i}) oc exp 



(3.6) 
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a homogeneous universe, one can assume this formation time to be independent of position, 
tf^i = tf, so that differences in ages between different redshifts can be used to estimate H{z) 
as in [41, 42, 44]. However, in an inhomogeneous cosmology, the formation time may depend 
on position within the void, which could significantly bias the resulting Hubble measurement. 
We discuss this issue in more detail in Appendix A (see also [46]). As a consequence, we need 
to treat the formation time at each position as an independent, free parameter^. 

Demanding only that the formation times tj^i > 0, and assuming a uniform prior, allows 
us to analytically marginalize over the formation times. 



where we again omit cosmology independent (additive) terms. 
3.4 The local matter density 

We use a measurement of the local matter density at redshift zero by [54], which we can treat, 
in the context of LTB models, as a measurement of Om(0) (equal to the parameter Ojn in 
the GBH parametrization and simply in ACDM). Using the mass-to-light ratio of galaxy 
clusters, they find Clm = 0.16 with uncertainty cjq^ = 0.06. We thus use the following 



More recent measurements of 0^ exist that give significantly smaller error bars (and 
values close to Qm = 0.28). However, these more recent determinations typically depend 
on the assumption of a homogeneous cosmology and often rely on a good understanding of 
the cosmic microwave background or the distribution of large scale structure. We choose the 
measurement by [54] because it depends largely on well understood local physics. It was 
derived by using the observed mass-to-light ratios of galaxy clusters to estimate the mass- 
to-light ratio of the local universe. This quantity was then multiplied by the observed total 
luminosity of the local universe to obtain the local matter density (we expect the cosmology 
dependence of the conversion factor between observed cluster mass-to-light ratio and total 
mass-to-light ratio, which was derived from simulations, to be small compared to the error 
bar). 

^Alternatively, if we could model the formation time's position dependence tf{r), the lost information on 
the Hubble rate can be restored, and information may even be added through potential cosmology dependence 
of tf{r). We will not pursue this approach here, as it will require modeling of structure formation, which we 
wish to avoid in this work. 




(3.7) 



This then gives us an effective of 






(3.9) 
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4 Results 



We now derive constraints from the data sets described in Section 3 (abbreviated as SN, Hq, 
Ages, Om), paying particular attention to model comparison between the void model and 
ACDM. For this purpose, we consider two statistics. First, we will compare the difference in 
best-fit values, corrected for the number of degrees of freedom, 

Xvoid - Xacdm - (dofvoid - dofACDM) = Ax^ - Adof, (4.1) 

with Adof = 4 — 2 = 2. Since neither model is nested inside the other, it is difficult to turn 
this statistic into a useful quantitative diagnostic to choose between the two models. 
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Figure 1. Marginalized posterior probability distributions of the void model parameters. 
Secondly, we will consider the Bayes factor, see e.g. [55, 56], 



B 



Eacd 



M 



where the Bayesian evidence is defined as 

E= I dep{e)p{D\e), 



(4.2) 



(4.3) 



and can be seen as the average likelihood of the data over a model's parameter space. Here, 
the integral is over the parameter space of the model, p{9) is a prior on the model parameters 
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6, and p(D\0) the likelihood given the data D. The Bayes factor can be interpreted in terms 
of betting odds when the models are deemed a priori equally likely, i.e. the odds are B : 1 
in favor of ACDM, and we will interpret it using the scale suggested by [55] (see also [57]), 
which we show in Table 1. Note that, while in the majority of cosmological applications the 
Bayes factor has been used to compare nested models, this statistic is of straightforward use 
and interpretation also for non-nested models, as is the case here. 



B 


evidence against void model 


1 to 3.2 


barely worth mentioning 


3.2 to 10 


substantial 


10 to 32 


strong 


32 to 100 


very strong 


> 100 


decisive 



Table 1. Jeffreys' scale for interpreting the Bayes factor. 

The Bayes factor strongly depends on the priors p{9). The (relevant part of) ACDM 
parameter space consists of and Hq and we will put the same uniform priors on these 
parameters as on their void model counterparts and Hqq. We choose the prior ranges on 
these parameters such that all significant likelihood is included within them, so that the ratio 
of evidences becomes independent of these ranges. The other void parameters R and AR do 
not have a ACDM counterpart and we choose priors R = — 50 Gpc and AR = 0.5 — 50 
Gpc, with an additional requirement that AR > O.IR for reasons of numerical stability (this 
latter requirement does not have a significant effect on the final probability distributions or 
Bayes factors). 

4.1 SN + Ho 

Considering first the supernova and local Hubble rate measurements, both models provide a 
good fit, as expected. Comparing the x^'s at the best-fit points in the two parameter spaces, 
we find Ax^ — Adof = 1.6, not strongly favoring either model. The Bayes factor, on the other 
hand, equals B = 13.2, thus (just) being in the range considered "strong" evidence for ACDM 
over the void cosmology, according to Jeffreys' scale. The reason that the Bayesian analysis 
more strongly disfavors the void model than the based analysis is that the Bayes factor 
penalizes models that need a large parameter space to obtain a good fit. The void model has 
two parameters more than ACDM, both {R and AR) being poorly constrained. 

The black curves in Fig. 1 show the marginalized probability distributions of the param- 
eters in the void model (see Fig. 2 for the ACDM case, also for the curves discussed in the 
following subsections). 

4.2 Adding galaxy ages 

In Fig. 3, we show the age of the universe as a function of redshift for the SN + Hq best- 
fit points in parameter space for both the void model and ACDM. Also plotted are the 
galaxy based red envelope ages described in Section 3.3, which provide a lower bound on 
the age of the universe. One sees that while the ACDM model (solid black line) is perfectly 
consistent with these observations, the void model (dashed black) is under some tension. More 
quantitatively, the difference in XAges between the two models (for these specific parameter 
values) is Axi = 48 (specifically, xiges = 49.9(2.3) for LTB (ACDM)). This difference in 



-9- 




Figure 2. Marginalized posterior probability distributions of the ACDM parameters. 

\^ is significantly alleviated once we allow the model parameters to vary (red curves) and we 
again compare best-fit values. The same tension can also be seen by comparing the black, 
SN + i^Oj curves in Fig. 1 to the red ones, which shows the parameter distribution preferred 
by the galaxy age data alone. Clearly, the galaxy ages prefer lower values of //qo (and f^in) 
than the SN + data, and require larger void dimensions. 

To get a complete picture of the effect of the galaxy age data, we next calculate the 
likelihoods for the full parameter spaces with the combined data set SN + + Ages. The 
resulting parameter constraints are shown by the blue curves of Fig. 1. 

Both models still provide a good fit to the data in the sense that the best-fit values 
are compatible with the number of observables minus number of degrees of freedom. However, 
the difference between the two models in best- fit x^'s is now A^^ — Adof = 5.4, showing that 
the void model is quite strongly disfavored relative to ACDM also based on the statistic. 
Comparing the Bayesian evidences gives B = 17, strengthening the evidence against the void 
model. 

4.3 Adding the local matter density 

As can be seen in Fig. 1, the combination of data discussed in the previous subsection requires 
a very low value of f^in in the void model. However, measurements of the local matter density 
exist that point at a larger value, more consistent with, the preferred. ACDM value of 
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Figure 3. Age of the universe versus redshift in void model and ACDM, compared to observed lower 
bounds from galaxy ages. When the void model parameters are chosen to provide an optimal fit to 
the supernova and local Hubble parameter measurements, the resulting cosmology is in significant 
conflict with the age observations (dashed black curve). ACDM, on the other hand, is consistent with 
the age data when the SN +Hq best-fit parameters are used (solid black). When the void model 
parameters are adjusted to agree with the age data and SN +Ho simultaneously, a much better fit to 
the age data can be reached (red dashed), but both the best-fit Ax^ statistic relative to ACDM and 
the Bayesian evidence change in favor of ACDM when these three data sets are combined (see text). 

Thus, using such measurements, the void model can be further disfavored relative to ACDM. 
As a proof of concept, we use the measurement described in Section 3.4, even though this 
measurement is rather old and its value is even low compared to ACDM (see Section 3.4 for 
a discussion). 

We show the parameter probability distributions preferred by this measurement alone 
in green in Fig. 1 (of course, it only constrains rijn and Qm) and the combined constraints 
from SN + Hq + Ages + Local Matter in magenta. The value of the best-fit Ax^ — Adof now 
equals 8.0 (although both models are still able to provide a good fit to the data in terms of 
the absolute values of the best fit x^'s) and the Bayes factor is B = 39.8, making the evidence 
in favor of a homogeneous universe "very strong" according to Jeffreys' scale. 

The best-fit void profiles for the different combinations of data sets are shown in Fig. 4. 
We show the final two dimensional 68% and 95% confidence level contours in Fig. 5. As 
can also be seen from Fig. 1, the parameters governing the void dimensions, R and AR, are 
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Figure 4. The best fitting density profiles to tlie tliree different data combinations considered in 
this section. The solid lines show the present physical matter density relative to the matter density 
outside the void as a function of radial position r. The dashed lines show the quantity ^mif) (see 
Eq. (2.10)), which gives the average matter density within r relative to a critical density defined in 
terms of the transverse Hubble parameter at r. While the supernova and Hq data allow for a void 
of "only" a few Gpc in radius, the galaxy age data prefer a larger void, several tens of Gpc in radius, 
although the size distribution is broad (see e.g. Fig. 5). 



clearly poorly constrained. Even with our broad priors, with upper limits of 50 Gpc, we 
do not enclose all the likelihood within the prior ranges. However, we do not consider even 
broader priors for computational reasons. We do note that if we broadened the priors, the 
evidence for the void model would decrease, favoring ACDM even more strongly than with 
the current prior. By tightening the prior ranges in R and Ai2, it is in principle possible to 
improve the Bayesian evidence for the void model and thus make the Bayes factor smaller. 
However, optimizing priors to get a more favorable result is against the spirit of Bayesian 
analysis, which is why we do not quote these numbers. We did check that we can at most get 
a 30% increase in the evidence, and thus about a 30% decrease in the Bayes factor. Therefore, 
optimizing the prior ranges would not cause a significant change in our main conclusions. 
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Figure 5. Marginalized two dimensional posterior probability distributions of the void model 
parameters. Black contours indicate 68% and 95% confidence level regions. The black dot indicates 
the best-fit point. The region below the red line in the lower part of the R — Ai? plane is excluded 
by our prior. 



5 Discussion and conclusions 

In this work, we have constrained LTB models proposed to explain the apparent accelerated 
expansion of the universe without introducing a cosmological constant. We have tested the 
LTB model, together with ACDM, against several types of data, including a recent compila- 
tion of supernova data, and the strongest constraint to date on the local expansion rate. New 
to this analysis is the addition of age data of old, passively evolving galaxies. We have argued 
that, while it is not allowed to use the differential ages as a Hubble parameter measurement 
in the context of void models, the data still have strong constraining power when considered 
as a set of lower bounds on the age of the universe in a wide redshift range z ~ — 1.9. 
Finally, we added a cluster-based measurement of the local (relative) matter density Vt^ to 
further constrain the models. 

In the previous section, we found that, while in terms of goodness of fit, both the void 
model and ACDM provide a satisfactory fit to all data combinations considered, a comparison 
of the relative merits of each model strongly favors the latter. Specifically, ACDM provides 
a better best-fit x^, with the difference (Ax^ = Xvoid ~ Xacdm)' ^-fter correcting for the fact 
that the void model has two more free parameters, being Ax^ + 2 = 1.6, 5.4 and 8.0 for SN 
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+ Hq, SN + i^o + Ages and SN + Hq + Ages + respectively. Bayesian analysis provides 
a good way of comparing the two models against the data through the Bayes factor, which 
is the ratio of Bayesian evidences. In addition to how well each model is able to fit the data, 
the advantage of this statistic is that it also takes into account Occam's Razor by penalizing 
models that need a large parameter space to obtain a good fit. The Bayes factors (here, the 
ratio of ACDM over void evidence) are B = 13.2, 17.0 and 39.8 for the data combinations 
listed above. Therefore, the first two combinations of data sets provide "strong" evidence 
against the void model according to Jeffreys' scale, while the combination of all data even 
gives "very strong" evidence. 

Our results are consistent with studies based on the kinetic SZ effect [24-28], which 
find that void models that can fit the supernova data are ruled out because they predict 
a much larger kSZ signal than consistent with observation (although this does assume the 
matter and radiation perturbations are adiabatic), and with studies combining the CMB 
temperature spectrum with other data sets. For instance, [31] shows, among other things, 
that the CMB in combination with supernovae require a very low value of Hqq, clearly in 
conflict with the measurement by [51], and that void models predict a very low clustering 
amplitude at low redshift, which is also inconsistent with data from clusters and other probes. 
The novel thing about our work is that it provides independent confirmation of the problems 
with large void models, without using CMB or LSS data. Moreover, our study suggests that 
combining the age data we used, with the data used in other studies, including the CMB and 
LSS data, could lead to significantly larger Bayes factors, even more firmly ruling out the void 
model. Of course, in that approach one would lose the independence of the details of cosmic 
perturbations in void models, which was one of the advantages of our approach in this work. 

Large values of the parameters governing the void size, R and AR, are not well con- 
strained by our data, with void sizes > 50 Gpc still allowed. These parameters have much 
stronger upper limits when CMB and LSS data are included. For instance, using SN + Hq + 
BAO + CMB data, [33] finds R and AR to be at most a few Gpc, with R = O.lStoit Gpc 
and AR = 2.56^j_l%^ Gpc at 68% confidence level. 

There are some potential caveats to the analysis presented here. First of all, we re- 
stricted ourselves to asymptotically flat cosmologies, Oout = 1- Li addition, we assumed a 
homogeneous big bang time such that at early times, the void is just a small perturbation 
onto a homogeneous Einstein-de Sitter cosmology. Finally, we did not consider the most gen- 
eral void profiles, instead relying on a specific parametrization involving four free parameters. 
Specifically, we did not include void profiles of the "unconstrained" type, which (partially) 
compensate the central underdensity with an overdensity near the edge of the void. 

Our study thus strongly constrains the class of the simplest, perhaps most reasonable 
LTB models. Relaxing some or all of the above assumptions might improve how well the 
model fits the data, but also adds more parameter space. It is thus not clear if this would 
strengthen or weaken the evidence against LTB cosmologies. Our main goal here was to 
introduce an alternative, strong method to constrain void models, and a study of constraints 
on the most general void cosmologies is beyond the scope of this article. 

In conclusion, void models as an explanation for the apparent cosmic acceleration are 
under attack from a number of independent types of cosmological data and are by now looking 
less and less like a viable alternative for a cosmological constant or dark energy. 
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A Inhomogeneous galaxy formation time 

The galaxy age data used in this work (see section 3.3) have in the past been used to estimate 
the Hubble parameter as a function of redshift. This method works because the passively 
evolving galaxy samples employed to determine the "red envelope" ages can in a homogeneous 
universe be assumed to have formation times whose statistics do not depend on position 
(i.e. on the redshift at which they are observed). The age difference between different redshifts 
then provides an unbiased estimator of the difference in cosmic time between those redshift, 
allowing a Hubble measurement through the relation 

H-' = - „ (A.l) 
dln{l + z) ^ ' 

In an inhomogeneous cosmology, this relation becomes 

SO that a differential cosmic time measurement could still carry valuable information. The 
problem with using ages is that, in an inhomogeneous universe, the formation time corre- 
sponding to red envelope ages at different redshifts may depend on position. 

How significant this position dependence is, is determined by the void profile at the time 
of galaxy formation. If at these early times the void is reduced to a small perturbation onto a 
homogeneous universe (which in the t — ?■ limit is the case) , any spatial variation of formation 
time can likely be neglected, but if the void still represents a significant underdensity, this 
variation may render the galaxy age-based Hubble measurement invalid. Since we do not 
know the typical formation times of the galaxies defining the red envelope, we consider in 
Fig. 6 a range of cosmic times^. The figure shows the void profile, quantified by the matter 
density relative to the matter density outside of the void, for the best-fit void profile to the 
data (SN + Hq + Ages + ^m), and shows that, even if typical formation times are as low 
as 0.3 Gyr, the underdensity is still approximately 30% in the void center. The later the 
formation time, the deeper the void is. 

We next study the effect on the Hubble parameter determination more quantitatively. 
Using 

t = tf + Age, (A.3) 

where t is cosmic time and tf the formation time, the usual estimator for the Hubble parameter 
would lead to a bias if tf = tf(r) because 



dln{l + z) ^ dln{l + z) 



To put these times in context, the average best-fit value of tf^i based on the full SN + Ho + Ages + flm 
data in the void model is about 1 Gyr, while in the ACDM case, we find 2.4 Gyr. A previous study of the 
effect of inhomogeneous formation time, [46], considered 0.3 Gyr as a plausible formation time. 
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Figure 6. The matter density relative to the matter density outside of the void at different potential 
galaxy formation times. The best-fit void parameters to the full set of data (SN + Hq + Ages + il,„) 
is used. Even at cosmic times as early as t — 300 Myr, the void underdensity is significant, with a 
~ 30% underdensity at the center. 



Thus, if tf{r) is known, one can estimate the relative bias in the (inverse) Hubble parameter. 



Hl 



dtfjrjz)) 
dln(l + z) 



(A.5) 



Since we do not have a good way of calculating tf(r), we consider a toy model, which, 
while not quantitatively correct, will at least tell us whether or not the bias is likely to be 
significant. The toy model simply assumes that tf{r) is reached when the local matter density 
reaches a critical value p^,, i.e. 

p{r,tf{r)) = p,. (A.6) 



This results in 



dt 



f 



dr 



d\n{l + z) 



dln(l + z) 



dp/dr 
dp/dt 



t=tfir) 



(A.7) 



which, together with Eq. (A.5), allows the calculation of AHJ^^ /H^^ once a void model is 
specified. We calculate this bias, assuming the best-fit void model (again, to SN + Hq + 
Ages + r^rn.)) for several choices of the formation time at the center of the void, and show the 
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results in Fig. 7. The bias in the toy model clearly depends strongly on the central formation 
time, with a formation time of 0.3 Gyr leading to a negligible bias compared to current error 
bars, but a formation time of 3 Gyr causing a bias of up to ~ 25%. 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Z 



Figure 7. The error in the determination of the (inverse of the) longitudinal Hubble parameter due 
to assuming the statistics of the galaxy formation time to be independent of location in the void. The 
calculation uses the best-fit void model to the SN + iJg + Ages + data set, and a simple model 
that assumes the galaxy formation time i/ is reached when some critical energy density is reached 
(see text for details). Depending on the time of formation in the center of the void, the error in iJ^ 
can be large, up to 15% at low redshift if the formation time 'w,tf ~ \ Gyr. 

Taking into account the large uncertainty in the central formation time, and in the 
theoretical evaluation of the spatial variation tf(r~) (our estimate being based on a mere toy 
model), we conclude that the bias in the Hubble parameter may well be significant and it is 
not safe to use the method by [43] without taking the effect of inhomogeneous formation time 
into account. For this reason, we have in this work only used the red envelope ages as lower 
bounds on the age of the universe, without estimating Ki^. If a better understanding of the 
formation of galaxies in void cosmologies were to lead to accurate predictions of t/(r) (perhaps 
with a small number of free parameters), the information on Hi^ could be partially or fully 
restored. In fact, in our toy model, the effect of tf{r^ on the differential age measurement is 
(void) cosmology dependent so that the spatial variation of formation time would potentially 
even add new information. 
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